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Abstract. Merging the Chandra and XMM-Newton deep surveys with the previously 
identified i? 05^ T surveys a unique sample of almost 1000 AGN-1 covering five orders 
of magnitude in 0.5-2 keV flux limit and six orders of magnitude in survey solid angle 
with ~ 95% completeness has been constructed. The luminosity-redshift diagram is 
almost homogeneously filled. AGN-1 are by far the largest contributors to the soft 
X-ray selected samples. Their evolution is responsible for the break in the total 0.5-2 
keV source counts. The soft X-ray AGN-1 luminosity function shows a clear change 
of shape as a function of redshift, confirming earlier reports of luminosity-dependent 
density evolution for optical quasars and X-ray AGN. The space density evolution 
with redshift changes significantly for different luminosity classes, showing a strong 
positive evolution, i.e. a density increase at low redshifts up to a certain redshift and 
then a flattening. The redshift, at which the evolution peaks, changes considerably 
with X-ray luminosity, from z KiO.5-0.7 for luminosities log Li,=42-43 erg s~^ to z « 2 
for log I/a; =45-46 erg s~^. The amount of density evolution from redshift zero to the 
maximum space density also depends strongly on X-ray luminosity, more than a factor 
of 100 at high luminosities, but less than a factor of 10 for low X-ray luminosities. For 
the first time, a significant decline of the space density of X-ray selected AGN towards 
high redshift has been detected in the range log La; =42-45 erg s~^, while at higher 
luminosities the survey volume at high-redshift is still too small to obtain meaningful 
densities. A comparison between X-ray and optical properties shows now significant 
evolution of the X-ray to optical spectral index for AGN-1. The constraints from 
the AGN luminosity function and evolution in comparison with the mass function of 
massive dark remnants in local galaxies indicates, that the average supermassive black 
hole has built up its mass through efficient accretion (e ~ 10%) and is likely rapidly 
spinning. 



1 Introduction 



In recent years the bulk of the extragalactic X-ray background in the 0.1-10 
keV band has been resolved into discrete sources with the deepest ROSAT, 
Chandra and XMM-Newton observations ,32 55 26 27 33 2 703^ . Optical identi- 
fication programmes with Keck |61l40l5l7| and VLT [65122) find predominantly 
unobscured AGN-1 at X-ray fluxes Sx > 10""'^'' erg cm~^ s~^, and a mixture 
of unobscured and obscured AGN-2 at fluxes 10"^"* > Sx > 10^^^'^ erg cm~^ 
with ever fainter and redder optical counterparts, while at even lower X- 
ray fluxes a new population of star forming galaxies emerges |35I58I1I36I56I8| . 
At optical magnitudes R>24 these surveys suffer from large spectroscopic in- 
completeness, but deep optical/NIR photometry can improve the identification 
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completeness significantly, even for the faintest optical counterparts |75l43j . A 
recent review article jj^l summarizes the current status of X-ray deep surveys. 

The AGN/QSO luminosity function and its evolution with cosmic time are 
key observational quantities for understanding the origin of and accretion history 
onto supermassive black holes, which are now believed to occupy the centers of 
most galaxies. X-ray surveys are practically the most efficient means of finding 
active galactic nuclei (AGNs) over a wide range of luminosity and redshift. Enor- 
mous efforts have been made by several groups to follow up X-ray sources with 
major optical telescopes around the globe, so that now we have fairly complete 
samples of X-ray selected AGNs. The most complete and sensitive sample was 
compiled recently by Hasinger, Miyaji and Schmidt |34| . concentrating on unab- 
sorbed (type-1) AGN selected in the soft (0.5-2 keV) X-ray band, where due to 
the previous ROS AT woik [5015 1| complete samples exist, with sensitivity limits 
varying over five orders of magnitude in flux, and survey solid angles ranging 
from the whole high galactic latitude sky to the deepest pencil-beam fields. These 
samples allowed to construct luminosity functions over cosmological timescales, 
with an unprecedented accuracy and parameter space. 



Table 1. The soft X-ray sample 



Survey'' 


Sohd Angle 


>S'xi4,lim 


iVtot 


A'agn-i'^ 






[dcg2] 


[cgs] 








RBS 


20391 


250 


901 


203 





SA-N 


684.0-36.0 


47.4-13.0 


380 


134 


5 


NEPS 


80.7-1.78 


21.9-4.0 


262 


101 


9 


RIXOS 


19.5-15.0 


10.2-3.0 


340 


194 


14 


RMS 


0.74-0.32 


1.0-0.5 


124 


84 


7 


RDS/XMM 0.126-0.087 


0.38-0.13 


81 


48 


8 


CDF-S 


0.087-0.023 


0.022-0.0053 


293 


113 


1 


CDF-N 


0.048-0.0064 0.030-0.0046 


195 


67 


21 


Total 






2566 


944 


57 



Abbreviations - RBS: The BOS'^T Bright Survey [Sll; SA-N: i?©^^!^ Selected Areas 
North [H; NEPS: ROSAT North Ecliptic Pole Survey [211; RIXOS: ROSAT Interna- 
tional X-ray Optical Survey g7|, RMS: ROSAT Medium Deep Survey, consisting of 
deep PSPC pointings at the North Echptic Pole P, the UK Deep Survey AS', the 
Marano field El and the outer parts of the Lockman Hole |6il39l : RDS/XMM: ROSAT 
Deep Survey in the central part of the Lockman Hole, observed with XMM-Newton 
|40l4ij|l<)| : CDF-S: The Chandra Deep Field South I65I75I43I : CDF-N: The Chandra 
Deep Field North [SEl. 
^ Excluding AGNs with z < 0.015. 

Objects without redshifts, but hardness ratios consistent with type-1 AGN. 
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2 The X-ray selected AGN— 1 sample 

For the derivation of the X-ray luminosity function and cosmological evolution 
of AGN well-defined flux-limited samples of active galactic nuclei have been 
chosen, with flux limits and survey solid angles ranging over five and six orders 
of magnitude, respectively (see Table 1). To be able to utilize the massive amount 
of optical identification work performed previously on a large number of shallow 
to deep iZ 05^47" surveys, the analysis was restricted to samples selected in the 
0.5-2 keV band. In addition to the i? 05^ T surveys already used in j5()l51| . data 
from the recently published ROSATNovth Ecliptic Pole Survey (NEPS) |29l54j . 
from an XMM-Newton observation of the Lockman Hole 021 as well as the 
Chandra Deep Fields South (CDF-S) |f)hl'7HI4:^| and North (CDF-N) |S|7| were 
included. In order to avoid systematic uncertainties introduced by the varying 
and a priori unknown AGN absorption column densities only unabsorbed (type- 
1) AGN, classified by optical and/or X-ray methods were selected. We are using 
here a definition of type-1 AGN, which is largely based on the presence of broad 
Balmer emission lines and small Balmer decrement in the optical spectrum of the 
source (optical type-1 AGN, e.g. the ID classes a, b, and partly c in [HI], which 
largely overlaps the class of X-ray type-1 AGN defined by their X-ray luminosity 
and unabsorbed X-ray spectrum j65| . However, as Szokoly et al show, at low X- 
ray luminosities and intermediate redshifts the optical AGN classification often 
breaks down because of the dilution of the AGN excess light by the stars in 
the host galaxy (see e.g. 52 ), so that only an X-ray classification scheme can 
be utilized. Schmidt et al. 61 have already introduced the X-ray luminosity in 
their classification. For the deep XMM-Newton and Chandra surveys in addition 
the X-ray hardness ratio was used to discriminate between X-ray type-1 and 
type-2 AGN, following 

Most (« 70-100%) of the extragalactic X-ray sources found in both the deep 
and wider X-ray surveys with Chandra and XMM-Newton are AGN of some 
type. Starburst and normal galaxies make increasing fractional contributions at 
the faintest X-ray flux levels, but even in the CDF-N they represent 20-30% 
of all sources (and create ^ 5% of the XRB). The observed AGN sky density in 
the deepest X-ray surveys is ~ 7200 deg"^, about an order of magnitude higher 
than that found at any other wavelength This exceptional effectiveness at 
finding AGN arises because X-ray selection (1) has reduced absorption bias 
and minimal dilution by host-galaxy starlight, and (2) allows concentration of 
intensive optical spectroscopic follow-up upon high-probability AGN with faint 
optical counterparts (i.e., it is possible to probe further down the luminosity 
function) . 

3 Number Counts and Resolved Background Fraction 

Based on deep surveys with Chandra and XMM-Newton, the X-ray log(N)- 
log(S) relation has now been determined down to fluxes of 2.4 x 10~^^ erg cm~^ 
s-\ 2.1 X lO"^*' erg cm-^ g-i, and 1.2 x IQ-^^ erg cm'^ s'^ in the 0.5-2, 2-10 
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Fig. 1. (a) Cumulative number counts N(>S) for the total sample (upper blue thin 
line), the AGN-1 subsample (lower black thick line), the AGN-2 subsample (red dotted 
line) and the galaxy subsample (green dashed line), (b) Differential number counts of 
the total sample of X-ray sources (open squares) and the AGN-1 subsample (filled 
squares). The dot-dashed lines refer to broken powerlaw fits to the differential source 
counts (see text). The dashed red line shows the prediction for type-1 AGN (from |34p. 

and 5-10 keV band, respectively |10I33I58I53I8| . Figure^ shows the normalized 
cumulative source counts N{> Sxi4)S^f^. The total differential source counts, 
normalized to a Euclidean behaviour (dN/d5'xi4 x S'^^^) is shown with open 
symbols in Figure QJd. Euclidean source counts would correspond to horizontal 
lines in these graphs. For the total source counts, the well-known broken pow- 
erlaw behaviour is confirmed with high precision. A broken power law fitted 
to the differential source counts yields power law indices of at, = 2.34 ± 0.01 
and aj = 1.55 ± 0.04 for the bright and faint end, respectively, a break flux 
of Sxi4 ~ 0.65 ± 0.10 and a normalisation of dN/dS'xi4 = 103.5 ± 5.3 deg~^ 
at S'xi4 = 1.0 with a reduced x^ = 1.51. We see that the total source counts 
at bright fluxes, as determined by the i?05'^T' All-Sky Survey data, are signif- 
icantly flatter than Euclidean, consistent with the discussion in [SI]- Moretti 
et al. on the other hand, have derived a significantly steeper bright flux 
slope {ab ~ 2.8) from ROSAT HKl pointed observations. This discrepancy can 
probably be attributed to the selection bias against bright sources, when using 
pointed observations where the target area has to be excised. 

The ROSAT HRI Ultradeep Survey had already resolved 70-80% of the extra- 
galactic 0.5-2 keV XRB into discrete sources, the major uncertainty being in the 
absolute flux level of the XRB. The deep Chandra and XMM- Newton surveys 
have now increased the resolved fraction to 85-100% |53l70j . Above 2 keV the 
situation is complicated on one hand by the fact, that the HEAO-1 background 
spectrum |46j . used as a reference over many years, has a ^ 30% lower normal- 
ization than several earlier and later background measurements (see e.g. |53|). 
Recent determinations of the background spectrum with XMM-Newton and 
RXTE inZ] strengthen the consensus for a 30% higher normalization, indicating 
that the resolved fractions above 2 keV have to be scaled down correspondingly. 
On the other hand, the 2-10 keV band has a large sensitivity gradient across the 
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band. A more detailed investigation, dividing the recent 770 ksec XMM-Newton 
observation of the Lockman Hole into finer energy bins, comes to the conclusion, 
that the resolved fraction decreases substantially with energy, from over 90% 
below 2 keV to less than 50% above 5 keV [23 • 

Type-1 AGN are the most abundant population of soft X~ray sources. For 
the determination of the AGN-1 number counts we include those unidentified 
sources, which have hardness ratios consistent with AGN-1 (a contribution of 
~ 6%, see Table Figure ^ shows, that the break in the total source counts 
at intermediate fiuxes is produced by type-1 AGN, which are the dominant 
population there. Both at bright fluxes and at the faintest fluxes, type-1 AGN 
contribute about 30% of the X-ray source population. At bright fluxes, they have 
to share with clusters, stars and BL-Lac objects, at faint fluxes they compete 
with type-2 AGN and normal galaxies (see Fig. ^ and [S]). A broken power 
law fltted to the differential AGN-1 source counts yields power law indices of 
ab = 2.55 ± 0.02 and af — 1.15 ± 0.05 for the bright and faint end, respectively, 
a break flux of S'xi4 = 0.53 ± 0.05, consistent with that of the total source 
counts within errors, and a normalisation of of dN/dS'xi4 = 83.2 ± 5.5 deg~^ 
at Sxi4 = 1.0 with a reduced x^=1.26. The AGN-1 differential source counts, 
normalized to a Euclidean behaviour (dN/dS'xi4 x ^xm) is shown with filled 
symbols in Figure ^ 

4 The Soft X— ray Luminosity Function and Space 
Density Evolution 

Hasinger, Miyaji and Schmidt "M* have employed two different methods to derive 
the AGN-1 X-ray luminosity function and its evolution. The first method uses 
a variant of the method, which was developed in |50j . The binned luminosity 
function in a given redshift bin Zi is derived by dividing the observed number 
Nobs{Lx, Zi) by the volume appropriate to the redshift range and the survey 
X-ray fiux limits and solid angles. To evaluate the bias in this value caused 
by a gradient of the luminosity function across the bin, each of the luminosity 
functions is fitted by an analytical function. This function is then used to predict 
Nmdi{Lx, Zi). Correcting the luminosity function by the ratio Nobs/Nmdi takes 
care of the bias to first order. 

The second method uses unbinned data. Individual T^nax of the RBS sources 
are used to evaluate the zero- redshift luminosity function. This is free of the bias 
described above: using this luminosity function to derive the number of expected 
RBS sources matches the observed numbers precisely. In the subsequent deriva- 
tion of the evolution, i.e., the space density as a function of redshift, binning 
in luminosity and redshift is introduced to allow evaluation of the results. Bias 
at this stage is avoided by iterating the parameters of an analytical represen- 
tation of the space density function. Together with the zero-redshift luminsity 
function this is used to predict Nmod{Lx, Zi) for the surveys. The observed den- 
sities in the bins are derived by multiplying the space density value by the ratio 
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Fig. 2. The soft X-ray luminosity function of the type-1 AGN sample in different 
redshift shells for the nominal case as labelled. The error bars correspond to 68% 
Poisson errors of the number of AGNs in the bin. The best-fit two power-law model 
for the 0.015 < z < 0.2 shell are overplotted in the higher redshift panels for reference. 
The dotted and dashed lines give the best-fit PLE and LDDE models (from 34 ). 



Nobs{Lx,Zi)/Nmod{Lx,Zi). At this stage, none of the densities are derived by 
dividing a number by a volume. 

The other difference between the two methods is in the treatment of missing 
redshifts for optically faint objects. In the binned method, all AGN without 
redshift with R > 24.0 were assigned the central redshift of each redshift bin to 
derive an upper boundary to the luminosity function. In the unbinned method, 
the optical magnitudes of the RBS sources were used to derive the optical redshift 
limit corresponding to i? = 24.0. The Knax values for surveys (such as CDF-N) 
spectroscopically incomplete beyond R — 24.0 were based on the smaller of the 
X-ray and optical redshift limits. 
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Fig. 3. Comparison between the space densities derived with two different methods. 
The blue datapoints with error bars refer to the binned treatment using the Nobs/Nmdi 
method, the dashed horizontal lines corresponding to the maximum contribution of 
unidentified sources. The thin and thick red lines and dots refer to the unbinned method 
(from 123]). 



Figure 121 shows the luminosity function derived this way in different redshift 
shells. A change of shape of the luminosity function with redshift is clearly 
seen and can thus rule out simple density or luminosity evolution models. In 
a second step, instead of binning into redshift shells, the sample has been cut 
into different luminosity classes and the evolution of the space density with 
redshift was computed. Figure |31 shows a direct comparison between the binned 
and unbinned determinations of the space density, which agree very well within 
statistical errors. 

The fundamental result is, that the space density of lower-luminosity AGN-1 
peaks at significantly lower redshift than that of the higher-luminosity (QSO- 
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type) AGN. Also, the amount of evolution from redshift zero to the peak is much 
less for lower-luminosity AGN. The result is consistent with previous determi- 
nations based on less sensitive and/or complete data, but for the first time our 
analysis shows a high-redshift decline for all luminosities Lx < 10^^ erg s~^ 
(at higher luminosities the statistics is still inconclusive). Albeit the different 
approaches and the still existing uncertainties, it is very reassuring that the gen- 
eral properties and absolute values of the space density are very similar in the 
two different derivations in. 

A luminosity-dependent density evolution (LDDE) model has been fit to the 
data. Even though the sample is limited to soft X-ray-selected type-1 AGN, 
the parameter values of the overall LDDE model are surprisingly close to those 
obtained by Ueda et al. 2003 for the intrinsic (de-absorbed) luminosity function 
of hard X-ray selected obscured and unobscured AGN, except for the normal- 
ization, where Ueda et al. reported a value about five times higher. 

These new results paint a dramatically different evolutionary picture for low- 
luminosity AGN compared to the high-luminosity QSOs. While the rare, high- 
luminosity objects can form and feed very efficiently rather early in the Universe, 
with their space density declining more than two orders of magnitude at redshifts 
below z=2, the bulk of the AGN has to wait much longer to grow with a decline 
of space density by less than a factor of 10 below a redshift of one. The late 
evolution of the low-luminosity Seyfert population is very similar to that which 
is required to fit the Mid-infrared source counts and background and also 
the bulk of the star formation in the Universe j^I], while the rapid evolution of 
powerful QSOs traces more the merging history of spheroid formation 23 . 

This kind of anti-hierarchical Black Hole growth scenario is not predicted in 
most of the semi-analytic models based on Cold Dark Matter structure forma- 
tion models (e.g. |37I71 |'). This could indicate two modes of accretion and black 
hole growth with radically different accretion efficiency (see e.g. ^j). A self- 
consistent model of the black hole growth which can simultaneously explain the 
anti-hierarchical X-ray space density evolution and the local black hole mass 
function derived from the Mbh — cr relation assuming two radically different 
modes of accretion has recently been presented in j^Hj ■ 

5 Optical versus X-ray selection of AGN— 1 

The space density of soft X-ray selected QSOs from the Hasinger et al. sam- 
ple is compared to the one of optically-selected QSOs at the most luminous 
end in Fig. 0] The z < 2 number density curve for optically selected QSOs 
(Mfij < —26.0) is from the combination of the 2dF and 6dF QSO redshift surveys 
| 14| . The z > 2.7 number densities from [6D| and (20j have been originally given 
for Ho=50 km s~^ Mpc~^, fim = 1, = 0. Their data points have been con- 
verted to Ho=70 km s"^ Mpc~\ = 0.3, f2A = 0.7 and the Mb threshold has 
been re-calculated with an assumed spectral index of ao = —0.79 (/^ oc v"'"), fol- 
lowing e.g. inni- The plotted curve from |50I20) is for Mb < —26.47 under these 
assumptions. A small correction of densities due to the cosmology conversion 
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Fig. 4. Comparison of tlie space density of luminous QSOs between optically selected 
and soft X-ray selected samples (from ,34j). The X-ray number densities are plotted 
for the luminosity class log Lx = 44 — 45, both for the binned and unbinned analysis 
with the same symbols as in Fig. El The dashed lines represent the one sigma range for 
Mj,^ . < —26.0 from [141 . multiplied by a factor of 16 to match the X-ray space density 
at z=2. The triangles at z > 2.7 with la errors are from '60' (SSG95) and !20j after 
a cosmology conversion (see text) and a scaling by a factor of 40 to match with the 
soft X-ray density at z ~ 2.7. As discussed in this paper, both the rise and the decline 
of the space density, behavior changes with Lx and therefore that the comparison can 
only be illustrative. 

causing redshift-dependent luminosity thresholds has also been made, assuming 
d(l>/dlogLB oc Lg^'^ |2U^ . The space density for the soft X-ray QSOs for the 
luminosity class 44 < logLx < 45 has been plotted both for the binned and 
unbinned determination. The Groom et al. ^] space density had to be scaled 
up by a factor of 16 in order to match the X-ray density at z ~ 2. The Schmidt, 
Schneider & Gunn / Fan et al. data points have been scaled by a factor of 40 
to match the soft X-ray data at z — 2.7 in the plot. There is relatively little 
difference in the density functions between the X-ray and optical QSO samples, 
although we have to keep in mind, that both the rise and the decline of the 
space density is varying with X-ray luminosity, so that this comparison can only 
be illustrative until larger samples of high-redshift X-ray selected QSOs will be 
available. 



10 Giinther Hasinger 




Fig. 5. (a): comparison of X-ray fluxes and AB2500 UV magnitudes for tiie sample of 
~ 1000 X-ray selected type-1 AGN from [SH (blue points) witli the ROSAT-observed 
optically selected SDSS QSOs [5^. Filled red squares give the standard SDSS QSO, 
while open red squares give the speciflcally selected high-redshift SDSS sample (see 
|69|). (b): Monochromatic 2 keV X-ray versus 2500 A UV luminosity for the same 
samples. The blue (dark) solid line shows a linear relation between the two luminosities, 
while the yellow (light) solid line gives the non-linear relation Lx ^ Lify" from the 
literature |69| . 

As a next step we directly study the X-ray and optical fluxes and luminosities 
of our sample objects and compare this with the optically selected QSO sample 
of Vignali et al. based on SDSS-selected AGN serendipitously observed in 
ROSAT PSPC pointings. Because of the inhomogeneous nature and different 
systematics of the different surveys entering our sample, the optical/UV mag- 
nitudes of our objects have unfortunately much larger random and systematic 
errors and are based on fewer colours than the excellent SDSS photometry. In 
our preliminary analysis we therefore calculated the AB2500 magnitudes simply 
extrapolating or interpolating the observed magnitudes in the optical filters clos- 
est to the redshifted 2500 A band, assuming an optical continuum with a power 
law index of -0.7, i.e. not utilizing the more complicated QSO spectral templates 
including emission lines which have been used in 69 . A spectroscopic correction 
for the host galaxy contamination, as done for the SDSS sample, was also not 
possible for our sample, however, for a flux and redshift-selected sample of 94 
RBS Seyferts j59) we have morphological determinations of the nuclear versus 
host magnitudes (see below). In all other aspects of the analysis we follow the 
Vignali et al. treatment. FigurelSJi shows 0.5-2 keV X-ray fluxes versus AB2500 
magnitudes for our sample objects (blue stars) in comparison with the Vignali et 
al. SDSS sample (X-ray detections are shown as filled red squares, upper limits 
as down-pointing triangles). It is obvious, that our multi-cone survey sample 
covers a much wider range in X-ray and optical fluxes than the magnitude- 
limited SDSS sample. Unlike the SDSS sample, our sample shows a very clear 
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correlation between X-ray and optical fluxes, but also a wider scatter in this 
correlation. 

Figure shows the monochromatic X-ray versus UV luminosity for the 
same data. Now the X-ray and optically selected samples cover a similar pa- 
rameter range at the high luminosity end, but the X-ray selected data reach 
significantly lower X-ray and UV luminosities than the optically selected sam- 
ple. Again, there is a larger scatter in the X-ray selected sample. The figure also 
shows two analytic relations between X-ray and UV luminosity: a linear rela- 
tion Lx ocLf/v and the non-linear behaviour Lx ocL^y found in the literature 
(e.g. ISni)- While the Vignali et al. optically selected sample clearly prefers the 
non- linear dependence (see also ^J), this is not true for the X-ray selected sam- 
ple, which is consistent with a linear relation, apart from the behaviour at low 
luminosities, where significant contamination from the host galaxy is expected. 




Redshift Redshift 

Fig. 6. X-ray to optical spectral index aox as a function of redshift for different lumi- 
nosity classes for the Hasinger et al. sample of ~ 1000 soft X-ray selected type-1 AGN. 
Blue (dark) stars give the values derived using the total integrated optical light, while 
the green triangles give the values using only the nuclear component from 93 RBS 
AGN, derived by Salvato |59| from detailed imaging decomposition. The solid black 
squares with error bars give the median and variance of the X-ray selected sample. 
Open red squares with error bars show the average aox values for the Vignali et al. 
optically selected QSOs. 



To check on any evolution of the optical to X-ray spectral index with redshift 
we calculated aox values, following 69^ for all our sample objects. In order to see 
possible luminosity-dependent evolution effectssimilar to those observed in the 
space density evolution, we divided our sample into the same luminosity classes 
as in Section ^ Figure shows the aox values determined for objects in the 
luminosity class 43-44 and 44-45, respectively, as a function of redshift. Apart 
from a few wiggles, which are likely due to the omission of the emission lines 
in the optical AGN continuum, there is no significant evolution with redshift. 
The optically selected sample, on the other hand, shows a significant trend with 
redshift and average values inconsistent with the X-ray selected sample for most 
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of the redshift range. The diagram also shows, that this discrepancy is hkely 
not caused by the missing host galaxy contamination correction in our analysis. 
From the relatively small number of nearby (z<0.1) of RBS sources, where a 
morphological fitting procedure has been used to subtract the host emission 
from the total magnitude |^ we can estimate the host dilution effect, which 
is clearly larger at lower X-ray luminosities and makes the discrepancy even 
larger. The immediate conclusion is, that the average optical to X-ray sample 
properties are dependent on systematic sample selection effects. 

6 X— ray Constraints on the Growth of SMBH 

The AGN luminosity function can be used to determine the masses of rem- 
nant black holes in galactic centers, using Soltan's continuity equation argument 
| 63| and assuming a mass-to-energy conversion efficiency e. For a non-rotating 
Schwarzschild BH, e is expected to be 0.054, while for a maximally rotating 
Kerr BH, e can be as high as 0.37 [Ml- The AGN demography predicted, that 
most normal galaxies contain supermassive black holes (BH) in their centers, 
which is now widely accepted (e.g. [211 and references therein). Recent deter- 
minations of the accreted mass from the optical QSO luminosity function are 
around 2eg"![ • 10^ MqMpc~^ |13l73j . Estimates from the X-ray background spec- 
trum, including obscured accretion power obtain even larger values: 6-9 ^Hl or 
8-17 P3 in the above units, and values derived from the infrared band [^01 or 
multiwavelength observations [HI are similarly high (8-9). Probably the most re- 
liable recent determination comes from an integration of the X-ray luminosity 
function. Using the Ueda et al. [HS] hard X-ray luminosity function including 
a correction for Compton-thick AGN normalized to the X-ray background, as 
well as an updated bolometric correction ignoring the IR dust emission, Marconi 
et al gSI derived paccr ~ S.Se^^ • IO^MqAIpc'^. 

The BH masses measured in local galaxies are tightly correlated to the galac- 
tic velocity dispersion |21I25| , and less tightly to the mass and luminosity of the 
host galaxy bulge (however, see 01]). Using these correlations and galaxy lu- 
minosity (or velocity) functions, the total remnant black hole mass density in 
galactic bulges can be estimated. Scaled to the same assumption for the Hubble 
constant (Ho=70 km Mpc~^), recent papers arrive at different values, mainly 
depending on assumptions about the intrinsic scatter in the BH-galaxy correla- 
tions^// = (2.4±0.8), (2.9±0.5) and {A.Gtlih^Q- 10^ MqMpc'^, respectively 
|,'-{l73l45j . The local dark remnant mass function is fully consistent with the above 
accreted mass function, if black holes accrete with an average energy conversion 
efficiency of e = 0.1 0^1 > which is the classically assumed value and lies between 
the Schwarzschild and the extreme Kerr solution. However, taking also into ac- 
count the widespread evidence for a significant kinetic AGN luminosity in the 
form of jets and winds, it is predicted, that the average supermassive black hole 
should be rapidly spinning fast (see also jl7l73| 1. Recently, using XMM-Newton, 
a strong relativistic Fe Ka line has been discovered in the average rest-frame 



Clues Irom Deep X-ray Surveys 



13 



spectra of AGN-1 and AGN-2 jM], which can be best fit by a rotating Kerr 
solution consistent with this conjecture. 
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Miyaji and Niel Brandt for the excellent cooperation on studies of the X-ray 
background. 

References 

1. Alexander D.M., Aussel, H., Bauer F.E., et al., ApJ, 568, L85 (2002) 

2. Alexander D.M., Bauer F.E., Brandt W.N., et al., AJ, 126, 539 (2003) 

3. AUer, M. C, Richstone, D., AJ 124, 3035 (2002) 

4. Appenzeller, L, Thiering, I., Zickgraf, F.-J., et al., A&A (Suppl.) 364, 443 (1998) 

5. Barger A. J., Cowie L. L., Mushotzky R. F., Richards, E. A., AJ 121, 662 (2001a) 

6. Barger A J, Cowie L L, Bautz M W, Brandt W N, Garmire G P, et al., AJ, 122, 
2177 (2001b) 

7. Barger A. J., Cowie L.L., Capak R, et al., AJ, 126, 632 (2003) 

8. Bauer F.E., Alexander D.M., Brandt W.N., et al., AJ in press, |astro-ph/040800T1 

(2004) 

9. Bower R. C, Hasinger, C, Castander, F. J. et al., MNRAS 281, 59 (1996) 

10. Brandt W. N., Alexander, D. M., Hornschemeier, A. E., Garmire, G. P., Schneider, 
D. P., et al., AJ 122, 2810 (2001) 

11. Brandt, W.N., Vignali C, Lehmer B.D., et al, this volume (2005) 

12. Brandt, W. N., Hasinger, G., ARAA in press (2005) 

13. Chokshi, A., Turner, E. L., MNRAS 259, 421 (1992) 

14. Groom, S. M., Smith, R. J., Boyle, B. J., Shanks, T., Miller, L., et al., MNRAS 
349, 1397 (2004) 

15. De Luca, A., Molendi, S., A&A 419, 837 

16. Duschl, W.J., Strittmatter, P.A., in Active Galactic Nuclei: from Central Engine 
to Host Galaxy Abstract Book, meeting held in Meudon, France, July 23-27, 2002, 
Eds.: S. Collin, F. Combes and 1. Shlosman. To be published in ASP Conference 
Series, p. 76 (2002) 

17. Elvis, M., Risaliti, G., & Zamorani, G., ApJ 565, L75 (2002) 

18. Fabian, A. C., Iwasawa, K., MNRAS 303, 34 (1999) 

19. Fadda, D., Flores, H., Hasinger, G., et al., A&A 383, 838 (2002) 

20. Fan, X., Strauss, M. A., Schneider, D. P., Gunn, J. E., Lupton, R. H., et al., AJ 
121, 54 (2001) 

21. Ferrarese, L., Merritt, D., ApJ 539, L9 (2000) 

22. Fiore F., Brusa M., Cocchia F., et al., A&A 409, 79 (2003) 

23. Franceschini, A., Hasinger, G., Miyaji, T., Malquori, D., MNRAS 310, L15 (1999) 

24. Franceschini A., Braito V. & Fadda D., MNRAS, 335, L51 (2002) 

25. Gebhardt, K., Bender, R., Bower, G., Dressier, A., Faber S. M., et al., ApJ 539, 
L13 (2000) 

26. Giacconi R., Rosati P., Tozzi P., et al., ApJ, 551, 624 (2001) 

27. Giacconi, R., Zirm, A., Wang, J.X., et al., ApJS, 139, 369 (2002) 

28. Gilh, R., Cimatti, A., Daddi, E., et al, ApJ 592, 721 (2003) 

29. Gioia, I., Henry, J. P., MuUis, C.R., et al., ApJS 149, 29 (2003) 

30. Haehnelt, M. G., Kauffmann, G., in Black Holes m Binaries and Galactic Nuclei, 
ed. L. Kaper, E. P. J. van den Heuvel, & P. A. Woudt, p. 364 (Berlin: Springer) 
(2001) 



14 Gxinther Hasinger 



31. Hasinger G., Burg R., Giacconi R., et al., A&A 275, 1 (1993) 

32. Hasinger G., Burg R., Giacconi R., et al., A&A 329, 482 (1998) 

33. Hasinger G., Altieri B., Arnaud M., et al., A&A 365, 45 (2001) 

34. Hasinger, G., Miyaji, T., Schmidt, M., A&A submitted (2005) 

35. Hornschemeier A.E., Brandt W.N., Garmire G.R, et al, ApJ 541, 49 (2000) 

36. Hornschemeier A.E., Bauer F.E., Alexander D.M., et al., AJ 126, 575 (2003) 

37. Kauffmann, G. & Haehnelt, M., MNRAS 311, 576 (2000) 

38. Kormendy, J., Gebhardt, K., AIP conference proceedings 586, 363 (2001) 

39. Lehmann I., Hasinger G., Schmidt M., et al., A&A 354, 35 (2000) 

40. Lehmann I., Hasinger G., Schmidt M., et al., A&A 371, 833 (2001) 

41. Madau, P., Pozzetti, L., Dickinson, M., ApJ 498, 106 (1998) 

42. Mainieri V., Bergeron J., Rosati P., et al., A&A 393, 425 (2002) 

43. Mainieri V., Rosati P., Tozzi, P., et al, A&A submitted (2004) 

44. Marconi, A., Hunt, L. K. ApJ 589, L21 (2003) 

45. Marconi, A., Risaliti, G., Gilli, R., Hunt, L. K., Maiolino, R., Salvati, M., MNRAS 
351, 169 (2004) 

46. Marshall, F. E., Boldt, E. A., Holt, S. S., et al., ApJ 235, 4 (1980) 

47. Mason K. O., Carrera, F. J., Hasinger, G., et al., MNRAS 311, 456 (2000) 

48. McHardy I.M., Jones, L. R., Merrifield, M. R., et al., MNRAS, 295, 641 (1998) 

49. Merloni, A., MNRAS 353, 1035 (2004) 

50. Miyaji, T., Hasinger G., Schmidt M., A&A 353, 25 (2000) 

51. Miyaji, T., Hasinger, G., Schmidt, M., A&A 369, 49 (2001) 

52. Moran, E. C., Filippenko, A. V. & Chornock, R., ApJ 579, L71 (2002) 

53. Moretti, A., Campana, S., Lazzati, D., et al., ApJ 588, 696 (2003) 

54. Mullis, C. R., Henry, J. P., Gioia, 1. M., et al., ApJ submitted l |astro-ph/0408304t 
(2004) 

55. Mushotzky, R. F., Cowie, L. L., Barger, A. J., Arnaud, K. A., Nature 404, 459 
(2000) 

56. Norman, C., Ptak, A., Hornschemeier, A., et al., ApJ 607, 721 (2004) 

57. Revnivtsev, M., Gilfanov, M., Sunyaev, R., et al., A&A 411, 329 (2003) 

58. Rosati P., Tozzi P., Giacconi R., et al., ApJ 566, 667 (2002) 

59. Salvato, M., Dissertation, Potsdam University (2002) 

60. Schmidt, M., Schneider, D. P., Gunn, J. E., AJ 110, 68 (1995) 

61. Schmidt, M., Hasinger, G., Gunn, J. E., et al., A&A 329, 495 (1998) 

62. Schwope A., Hasinger, G., Lehmann, I., et al., AN, 321, 1 (2000) 

63. Soltan, A., MNRAS 200, 115 (1982) 

64. Streblyanska, A., Hasinger, G., Finoguenov, A. et al., A&A in press (2005) 

65. Szokoly G., Bergeron J., Hasinger G., et al., ApJS in press ( |astro-ph/0312324t 
(2004) 

66. Thorne, K. S., ApJ 191, 507 (1974) 

67. Tozzi, P., Rosati, P., Nonino, M., et al., ApJ 562, 42 (2001) 

68. Ueda, Y., Akiyama, M., Ohta, K., Miyaji, T., ApJ 598, 886 (2003) 

69. Vignah C., Brandt W.N., Schneider D.P., AJ 125, 433 (2003) 

70. Worsley, M., Fabian, A.C., Mateos, S., et al., MNRAS 352, L28 (2004) 

71. Wyithe, J.S.B., Loeb, A., ApJ 595, 614 (2003) 

72. Yang, Y., Mushotzky, R. F., Barger, A. J., et al., ApJ 585, L85 (2003) 

73. Yu, Q., Tremaine, S., MNRAS 335, 965 (2002) 

74. Zamorani, G., et al. A&A, 346, 731 (1999) 

75. Zheng W., Mikles, V. J., Mainieri, V., et al., ApJ in press, astro-ph/0406482 (2004) 



